Orbital Wave and its Observation in Orbital Ordered Titanates and Vanadates 
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We present a theory of the collective orbital excitation termed orbital wave in perovskite titanates 
and vanadates with the triply degenerate t2 g orbitals. The dispersion relations of the orbital waves 
for the orbital ordered LaVOa, YVO3 and YTiC>3 are examined in the effective spin-orbital cou- 
pled Hamiltonians associated with the Jahn- Teller type couplings. We propose possible scattering 
processes for the Raman and inelastic neutron scatterings from the orbital wave and calculate the 
scattering spectra for titanates and vanadates. It is found that both the excitation spectra and 
the observation methods of the orbital wave are distinct qualitatively from those for the e 9 orbital 
ordered systems. 
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I. INTRODCUTION 

When degenerate electron orbitals are partially filled 
in correlated electron systems, this is recognized to be 
an internal degree of freedom belonging to an electron 
as well as spin and charge. This orbital degree of free- 
dom has recently attracted much attention especially in 
transition-metal (TM) oxides. 1,2 In particular, the or- 
bital ordering (00) and fluctuation play a key role in 
anisotropic electric, magnetic and optical properties in 
correlated oxides. 

The collective orbital excitation in the orbital ordered 
state is known to be orbital wave (OW), in analogous 
to the spin wave in magnetically ordered state, and 
its quantized object is termed orbiton. The theoreti- 
cal study of OW has started in the idealized spin-orbital 
coupled model where the continuous symmetry exists in 
the orbital space, 3 and progressed in the anisotropic or- 
bital model. 4 A realistic calculation has been done in 
LaMn03 with the doubly degenerate e g orbitals. 5 Re- 
cently new peak structures observed in the Raman spec- 
tra in LaMn03 are interpreted as scatterings from OW. 6 
Although there are some debates about origin of the 
newly found peaks, OW or multiphonons, ' the en- 
ergy and polarization dependences of the observed Ra- 
man spectra show a good agreement with the calculation 
based on the OW interpretation. 

The perovskite titanates i?Ti03 and vanadates -RVO3 
(i?: rare-earth ion) with the triply degenerate ti g or- 
bitals are another class of materials where OW is ex- 
pected. One of the well known orbital ordered mate- 
rials is YTi03 where the nominal valence of the tita- 
nium ion is 3+ and one electron occupies the triply de- 
generate t2 g orbitals. The orbital ordered state associ- 
ated with the Jahn- Teller ( JT) type lattice distortion has 
been confirmed by the resonant x-ray scatterings, NMR 
and the polarized neutron scattering experiments and so 
on. 9-12 These results almost coincide with each other: 
there are four different orbitals in a unit cell where the 
wave functions are given as -j^(d xy + d yz ), -^(d xy - d yz ), 

-^{d xy + d zx ) and -^{d xy - d zx ). 13 This type of 00, 



termed the {d y(x+z) /d y(x _ z) /d x(y+z) /d x{y _ z) )-type from 
now on, is also supported by previous calculations. 14-16 
Although an exotic orbital state is proposed recently by 
taking into account the quantum fluctuation, 17 ' 18 the 
predicted orbital state being incompatible with the crys- 
tal lattice symmetry is different from the above type of 
00. 

A series of -RVO3 is systematically examined in the 
recent studies. 14 ' 15 ' 19-30 Two electron occupy the triply 
degenerate orbitals in a V 3+ ions. A sequential phase 
transition is found in YV0 3 ; 19-21 the G-type 00 (O-G) 
occurs at Tboi=200K and the C-type antiferromagnetic 
(AFM) ordering (S-C) appears at 7Vi=115K. With fur- 
ther decreasing temperature, another orbital and mag- 
netic transitions appear at TV 2 = 7oo2=77K where the 
C-type 00 associated with the G-type AFM order (the 
(S-G/O-C) order) is realized. On the other hand, in 
LaV0 3 , the C-type AFM ordering occurs at Tat(=143K) 
and, at slightly below this temperature, the G-type 00 
(the (S-C/O-G) order) appears. 20-22 Types of 00 in 
vanadates are determined that the d xy orbital is occupied 
at all the vanadium sites and the d yz and d zx orbitals are 
alternately ordered in the xy plane (the C-type 00), and 
in all direction (the G-type 00). 24 ' 25 These kinds of 00 
are termed the pure 00 states, and the 00 such as re- 
alized in YTi03 is termed the mixed 00 from now on. 
It is suggested that the experimentally observed type of 
00 in vanadates associated with the JT type distortion 
explains the several optical and magnetic properties. 28,29 

Here, we present a theory of OW in titanates and vana- 
dates where the orbital ordered states are confirmed ex- 
perimentally well. The dispersion relations of OW are 
examined in the (S-C/O-G)- and (S-G/0-C)-phases for 
LaV03, YVO3 (the low temperature phase), respectively, 
and the {d y{x+z) /d y{x _ z) /d x(y+z) /d x{y ^ z) )-type 00 for 
YTi03. The calculations are based on the effective spin- 
orbital coupled Hamiltonians associated with the JT type 
electron-lattice coupling. We propose possible scattering 
processes for the Raman and inelastic neutron scatterings 
from OW and calculate the spectra. It is found that the 
excitation spectra and the observation methods of OW 
arc distinct qualitatively from those for the e g orbital 
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ordered systems such as manganites. 

In Sec. II, the model Hamiltonian for titanates and 
vanadates with the triply degenerate t 2g orbitals are in- 
troduced. In Sec. Ill, the dispersion relations of OW for 
LaVC>3, YVO3 and YTiC>3 are examined. The scattering 
spectra for the Raman and inelastic neutron scatterings 
from OW are shown in Sec. IV. Section V is devoted to 
summary and discussion. 

II. MODEL HAMILTONIAN 

The collective orbital excitations in the orbital ordered 
state are studied in the spin-orbital model derived from 
the generalized Hubbard Hamiltonian with the triply de- 
generate t 2g orbitals: 

H = Ed^2, A~t<jdiiv + Ti-el-el 

+ E (ty'd\ ia d jlla +H. c ), (1) 

with the electron-electron interaction term 

Hel-el — ^E n *7T n *7i + U' X/ n n n n' 

+ I X d^ad^, a ,d ila i d^' a 

i,7>7' ,<7,<t' 

di ia is the annihilation operator for the t 2g electron at 
site i with spin a = (T,|) and orbital 7 = (yz, zx,xy). 
The number operators are defined by nj 7<7 = d\ di 1<7 
and rii 7 = ^^n^a-. U and U' are the intra- and inter- 
orbital Coulomb interactions, respectively, I is the ex- 
change interaction, and J is the pair-hopping interaction. 

in Eq. (1) is the electron transfer integral between 
site i with orbital 7 and nearest neighboring (NN) site j 
with 7'. In the ideal perovskitc structure, the hopping 

integral is simplified as tJJ = t 5 77 >(8 7= (ih) + <5 7 =(fc;)) 
where I indicates a direction of a bond connecting sites 
i and j, and (h,k,l) = (x,y,z), (y,z,x), (z,x,y). The 
electron hopping occurs through the O 2p orbitals in be- 
tween the NN TM sites. 

We derive the effective Hamiltonian for titanates where 
the nominal electron configuration of the TM ion is d 1 . 
This Hamiltonian is defined in the Hilbert space where 
the electron number per site is restricted to be one or 
less due to the strong on-site Coulomb interactions. The 
Hamiltonian is classified by the intermediate states of the 
perturbational processes as derived in Rcf. 13; 

Tin = + Ht 2 + He + , (3) 

with 




The superexchangc interactions are given by J r = 
tl/E^ ] (r = Ti, T 2 , E, Ai) where £^ 2) 's are the energies 
of the intermediate states: Ejf* = U' — I, E E ^ = U — I, 

E^ = U' + I and = U + 21 where the relations 
U = U' + 21 and I — J are used. Jt x is the largest 
among them. The spin degree of freedom is described 
by the operator S t = \ J2-y SS ' d\ ls a S s'di ls ' with the Pauli 
matrices a. The orbital degree of freedom is represented 
by the eight orbital operators 0^7 classified by the irre- 
ducible representations of the Oh group as (T7) = (Eu), 
(Ev), (T 2 x), (T 2 y), {T 2 z), (T lX ), {T iy ), (T lZ ). These 
operators are defined by the generators of the SU(3) Lie 
algebra, i.e. the 3x3 Gell-Mann matrices X m (I — m <~ 8) 
as 31 

^ ir 7 — ~?K X/ diucr^mapdipa, (8) 

CT,Q,/3 

where (r 7 ,m) = {Eu,8), (Ev,3) (T 2 x,6), (T 2 y,4), 
(T 2 z, 1), {T\x, 7), (Tiy, 5), (T±z, 2). The operators OiE 7 
and OiT 2 -y represent the diagonal and off-diagonal com- 
ponents of the electric quadrupole moments, respectively, 
and 0iTi7 does the magnetic dipole ones, i.e. the orbital 
angular momentum. By utilizing the above orbital oper- 
ators, the orbital parts of the Hamiltonian in Eqs. (4)- (7) 
are given as 

A l ± = Wlw]±O l iEv O l jEv , (9) 
B l = V*Wj + W-Vj, (10) 
C l ± = 2 {O lT2l O jT2 i ± O iTll O jTl i) , (11) 

with 

Wl = | - ^-0\ Eu , (12) 
Vl = \ + \j\o\ Eu . (13) 
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°\eu and °\ev are defined by 

( °\eu \ _ ( cos ir TO ' sin f m \ ( ° lEu \ n a) 
V °Iev )~ v - sin ¥ "»« cos J V °«*< / ' 

with to; = (1, 2, 3) for a direction of the bond I = (x, y, z). 
In analogy with the spin Hamiltonian, A l ± and B l cor- 
respond to the Si Z Sj Z term, and C l + and C l _ to the 
Si+Sj- + Si-Sj+ and Si+Sj+ + Si-Sj- terms respec- 
tively. The latter term originates form the pair hopping 
processes which break the conservation of the total elec- 
tron number at each orbital, e.g. ^(ni I9 ). 

The effective Hamiltonian for vanadates, where the 
nominal electron configuration is d , is derived in a sim- 
ilar way from Eq. (1). The d 2 state is assumed to be the 
lowest 3 Ti state due to the Hund rule. The explicit form 
is given by 

T~tv = 7~(-Ai + T~Le + Wji + ^t 2 ) (15) 

with 

Ua 1 = - Ja 1 ^2 t (2 + Si ■ Sj"j 

(ij) 

x (-A l + +A l _+B l - 2C l + ) , (16) 
He = -Je g (i ~ $i ' 3j) 

(ij) 

x (-A l + + A l _ + B l + C l + ) , (17) 
n Tx = -Jt x E \ (l - Si ■ Sj) (A C l _) , (18) 

(ij) 

n T2 - -Jt 2 ]T \ i 1 ^ ■ Sj) (A 1 + &_) . (19) 

(ij) 

We introduce the spin operator Si with a magnitude 
S = 1, and the exchange parameters defined by Jr = 
tl/AE r where AE r = - with AE Al =U'-I, 
AE E = AE Tl =U' + 2I and AEt 2 = U' + AI. It is noted 
that a similarity between Hti and Tiy is attributed to 
the fact that, with respect to the orbital degree of free- 
dom, the high spin 3 T\ state for the d 2 configuration 
in the hole picture is equivalent to the d 1 state in the 
electron picture. Similar types of the spin-orbital cou- 
pled Hamiltonian in the triply degenerate tig orbitals 
are also derived by several authors for the d 1 and d 2 
systems 16 ' 18 ' 27 ' 30 ' 32 ' 33 

In addition to the electronic Hamiltonian, the electron- 
lattice interactions are introduced. In the tig orbital sys- 
tems, there are two kinds of the JT type interactions; 

HjT = 9E E QiE-fOiE-y 
i^ — {u,v) 

+ 9t 2 Q^T 2l O lT2l , (20) 



where gs and gx 2 are the coupling constants and QiE-y 
and QiT 2 ~t are the normal modes in an Oq octahedron 
with symmetries E 9 and T2g, respectively Qej directly 
modifies the TM-0 bond lengths and Qr 2 -y modifies the 
O-TM-0 bond angles. 

Energy parameter values have been numerically eval- 
uated by several authors. The effective exchange param- 
eters J s in the Heisenberg model are obtain from the 
spin-wave dispersion relations as J,f=5.5meV (z axis) 
and J;f =5.8meV (xy plane) in the (S-G/O-C) phase 
in YVO3, Jf = 2.2 - 4mcV, J$ y =2.6mcV in the (S- 
C/O-G) phase in YV0 3 (Ref. 26), and jf = jJ;=3meV 
in YTiC>3 (Ref. 17). The effective exchange parameters 
for the orbital operators in LaVC>3 are also estimated 
as Jf ee 4tl/(U' - 7)=33meV and J°=2meV (Ref. 28). 
The JT stabilization energy is obtained from the LDA+U 
method 28 as EjT=27mcV which is comparable or larger 
than the exchange interactions. The relativistic spin- 
orbit interaction which is not taken into account in the 
present model is about 0.4meV being much smaller than 
both the exchange and JT energies. 15 This is consistent 
with the experimental results in the magnetic x-ray scat- 
tering in YTiOs; the angular momentum separately esti- 
mated from the spin momentum is found to be negligible 
small. 34 Thus, the Hamiltonian Wti(v) +~H.jt introduced 
above is the minimal model for examination of OW. 

Here we mention the implications of the theoretical 
model for the observed spin/orbital orders. As we have 
shown in Ref. 13, in the mean field theory, the large 
orbital degeneracy remains in the ferromagnetic (FM) 
ground state in Hti- A small perturbation, such as 
the Jahn- Teller type distortion, the relativistic spin-orbit 
coupling, the GdFeGvtype lattice distortion, lifts the de- 
generacy. We suppose that the Jahn- Teller type distor- 
tion with the Tig symmetry, gr 2 QiT 2 m Eq- (20), plays 
a key role to stabilize the observed mixed orbital order. 
The (S-C/O-G) order for LaV03 is reproduced by 7iy 
The FM (AFM) order along the c axis (in the ab plane) 
is attributed to the alternate (uniform) alignment of the 
d yz and d zx orbitals (the d xy orbital) along the c axis 
(in the ab plane). The (S-G/O-C) order for the YV0 3 
is obtained by Tiy and the Jahn- Teller type interaction 
with the Eg symmetry The AFM spin order and the uni- 
form alignment of the d yz {d zx ) orbital along the c axis is 
stabilized cooperatively, as discussed in the next section. 

III. ORBITAL WAVE 

The dispersion relations of OW in the orbital or- 
dered states are obtained by utilizing the Holstcin- 
Primakoff transformation for the generators in the SU(3) 
algebra. 35 ' 36 For example, at a site where the d xy orbital 
is occupied, there are two excitation modes; an excitation 
between the d xy and d yz orbitals denoted by a boson op- 
erator y (yt), and that between d xy and d zx denoted by 
x (2^). The orbital operators Oirj are transformed into 
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FIG. 1. Dispersion relations of OW in the (S-C/O-G) and 
(S-G/0-C)-phases for LaVOa and YVO3 (the low temper- 
ature phase), respectively. The absolute values of the spin 
correlation function between NN sites \{Si ■ Si+s) \ are 1 (bold 
lines), 0.5 (broken lines) and (dotted lines). The parame- 
ter values are chosen to be / /U'= 0.125. (?bQb/Jo=0.8 and 
gT 2 Qr 2 =0. The dispersion curves in the (S-G/O-C) phase 
and the curve in the (S-C/O-G) phase with \(Si ■ S i+ s} \ = 1 
are doubly degenerate. 



these boson operators as 
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1 

V2 



OiEu 
OiEv 

V 

OiT a z = ^ (y\xi ± x\y^j , 
O.T a y = Q (y/l - N m ± yjy/l-Ni) , 
O iT<xX = (~A (x\y/l-Ni± y/l - W iXi ) , 



(21) 



with a = (1,2) . The plus and minus signs in OiT aX , 
OiT aV and OiT aZ are for the a = 1 and a = 2 cases, re- 
spectively. We define Ni = n 



with rii X = x\xi 



and iiiy = y\y%. In the linear spin wave approxima- 
tion, \/l — Ni is replaced by 1. We have cheked that 
the Green's function method for the operators 0^7, i-e. 

Gv^'it-t^n-rv) = -ie(t-t)([Oiri(t),Oi, T >yn> 
with the decoupling approximation, reproduces the cal- 
culated results for OW with the boson method intro- 
duced above. 

In Fig. 1, we present the dispersion relations of OW 
in the (S-C/O-G) and (S-G/O-C) phases for LaV0 3 and 
YVO3 (the low temperature phase), respectively. The 
parameter values are chosen to be I/U' = i?=0.125, 
9eQe/Jo—0-S and gT 2 Qr 2 =0- The energy parameters 
are normalized by J = 4^/ (U' — I) which is estimated to 
be about 33meV for LaV0 3 (Ref. 28) and is supposed to 



be smaller in YT1O3 due to the larger GdFe03-type lat- 
tice distortion. The ratios of the exchange parameters are 
represented by the parameter R as Jt 2 / Jt ± — Je/ Jt x = 
(1 — J2)/(l + 2J2) and J Al /J Tl = (1 — i2)/(l + 4J2). Both 
the (S-C/O-G)- and (S-G/0-C)-phases, there arc four 
modes of OW attributed to the two different orbital oc- 
cupied sites in a unit cell; there are the excitations yA 
) and za (d zx — > d yz ) at site A where the 
d xy and d zx orbitals are occupied, and the excitations 
xb {d X y — > d zx ) and zb (d yz — > d zx ) at site B where the 
d xy and d yz orbitals are occupied. Two of the four, i.e. 
y A and xb, are the local modes within the linear spin 
wave theory and do not show dispersions. This charac- 
ter does not depend on the spin arrangements. These 
local modes originate from the facts that (1) the excited 
d xy hole does not hop along the z axis due to the orbital 
symmetry, and (2) a coherent motion of the excited d xy 
hole in the xy plane are impossible, since this motion is 
associated with increasing the number of the wrong or- 
bital arrangements. The latter implies that the orbital 
exchange processes do not recover the wrong orbital ar- 
rangements, and the triply degenerate orbital model is 
qualitatively different from the Heisenberg model with 
5=1. The remaining modes, z A and zb, are disper- 
sive along the z direction. The dispersion relation of 
OW in the (S-C/O-G) phase is explicitly obtained as 

E(k) = p Al tJ{K% + Kf)~ 2 ~ (pq cosak z ) 2 where we 
assume 1 = and qeQe — 9t 2 Qt 2 = 0. We introduce 
K x 2 [ z) = 2+(Si-Si +5x{6z) ) and Kf z) = l-(Si-S i+5x{5z) ). 
This energy has its minimum at k z = and the energy 
gap is attributed to the anisotropy in the orbital space, 
i.e. a lack of the SU(3) symmetry, in the orbital part 
of the Hamiltonian. In comparison with the OW in the 
(S-C/O-G) phase, the OW in the (S-G/O-C) phase is 
barely stable; with decreasing the spin correlation which 
corresponds to increasing temperature toward TV, a re- 
markable softening around k = (0, 0, it) occurs. In the 
case where 1 = and gT 2 Qr 2 — 0, the dispersion relation 

is given as E(k) = "^Ja^K^ — K*) cos ak z + \J\gEQE 
which has its minimum value at k z = ir, and the energy 
gap is attributed to the JT type interaction. This re- 
sult suggests an instability of the (S-G/O-C) phase to 
the (S-C/O-G) one with increasing temperature. This is 
consistent with the experimental fact that in i?V03 the 
O-C phase appears associated with the S-G order, and is 
changed into the (S-C/O-G) phase at 77K in YVO3. 21 

In Fig. 2, we present the dispersion relations of OW 
in the {d y{x+z) /d y{x _ z) /d x{y+z) /d x ( y _ z) )-type OO for 
YTi03. The FM and paramagnetic spin orders are as- 
sumed, and the exchange parameter is taken to be R = 
I/U' = 0.125. The JT type interaction parameters arc 
chosen to be qeQe/Jo = gT 2 Qr 2 /Jo = 1-2, where 9eQe 
is larger than that in vanadates and Qt 2 Qt 2 is introduced. 
This is based on a consideration that, in comparison with 
vanadates, J is supposed to be smaller due to the large 
GeFe0 3 -type lattice distortion, and Qt 2 is found to be 



4 




z r MX r A 



FIG. 2. Dispersion relations of OW 

in the (d y(x+z )/d y{x - z) /d x(y+z) /d x(y _ z) )-type OO for YTi0 3 . 
The bold and broken lines are for the FM and paramagnetic 
states, respectively. The parameter values are chosen to be 
I/U'=0.125, gsQsl 'Jo=9t 2 Qt 2 / 'J = 1.2. The Brillouin zone 
for the tetragonal symmetry is adopted. 

larger in the crystal structure of YT1O3. 12 > 22 > 23 In con- 
trast to the case of vanadates, all the eights modes, at- 
tributed to the four different orbitals in a unit cell, arc 
dispersive along all the directions in the Brillouin zone. 
This originates from the OO states with the mixed or- 
bitals where the excitations propagate along the three 
directions. The OW dispersions are classified into the 
two groups: for example, in the -j=(d xy + d yz ) orbital 
occupied site, the higher energy bands with smaller band 
widths are mainly attributed to the excitations to the d zx 
orbital, and the lower energy ones with a larger width are 
attributed to the excitations to the -j=(d xy — d yz ) orbital. 

IV. OBSERVATION OF ORBITAL WAVE 

A. raman scattering 

As introduced in Sec. I, in orbital ordered LaMn03, 
the new peak structures in the Raman spectra were ex- 
plained successfully as the scattering from OW. 6 Here 
we consider the Raman scattering as a possible probe 
to detect the OW in titanates and vanadates with the 
triply degenerate t 2g orbitals. It is considered the inter- 
site scattering processes where OW's are excited through 
the electronic exchange processes between the NN TM 
sites. 37 ' 38 This is attributed to the fact that the low- 
est electronic excitations in titanates and vanadates oc- 
cur across the Mott-Hubbard gap, unlike the mangan- 
ites where the electronic excitations across the charge- 
transfer gap dominate the excitation of OW. Depending 
on the types of OO, there are the following two scatter- 
ing processes: (i) The two-orbiton scattering: Consider a 
pair of the NN TM sites where the occupied orbitals are 
different and the pure orbitals (Fig. 3), such as the OO 
in vanadates. Through the second order processes of the 



(a) 




FIG. 3. Scattering processes in the Raman scattering, (a) 
and (b) The two-orbiton scattering processes in the pure OO 
states, (c) The one-orbiton scattering processes in the mixed 
OO state. 

interaction between photons and electrons, electrons at 
the two sites are exchanged and, at the final state, the oc- 
cupied orbitals are changed at both the sites. This is the 
analogous to the two magnon Raman scattering in the 
antifcrromagnets. Another two-orbiton process occurs in 
a pair of the NN TM site where the same kind of orbitals 
are occupied. In the intermediate state of the scattering 
process, where two electrons occupy the same orbital at 
a site, the occupied orbital is changed due to the pair- 
hopping interaction. At the final state, the two orbitons 
are created, (ii) The one-orbiton scattering: When the 
electrons occupy the so-called mixed orbital such as that 
in YTi03, the electron hops from one orbital (7) to the 
different orbital (7') in its NN TM site. When this elec- 
tron comes back to the orbital 7' in the initial site, one 
orbiton is excited. 

There are alternate two scattering processes from OW 
where one orbton is created at a TM site, (i) An electron 
is excited from a orbital 7 to 7' at the same site associ- 
ated with a creation of odd-parity phonons. Then, these 
phonons are annihilated by emitting a photon. Such 
kinds of the Frank-Condon processes have been consid- 
ered for the orbiton+phonon excitation in the optical 
conductivity spectra, and for the multiphonon excitation 
in the Raman spectra. 39 The total scattering-cross sec- 
tion ratio of this process to the inter-site process is of the 
order of 1CP 3 <~ 10~ 2 . (ii) The incident photon excites 
an electron from the 3<i 7 orbital to one of the Ap or- 
bitals at the same site. Then, this electron is relaxed to 
the 3d 7' orbital by emitting a photon. This scattering- 
cross section ratio to the inter-site process is estimated 
for titanates and vanadates to be of the order of 10 _1 . 

Here, we calculate the Raman spectra for the inter- 
site scattering process. It is supposed that this process 
provides the main contribution for the OW scattering, 
in particular, in the two-orbiton energy regions. The en- 
ergy, momentum and polarization of initial (scattered) 
photon are fcj(fc/) and Xi(Xf), respectively. The 

differential scattering cross section from OW is given as 
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FIG. 4. Raman spectra S XiX f (u) for the (S-C/O-G)- and 
(S-G/O-C) phase for LaVOs and YVO3 (the low temperature 
phase), respectively. The two-orbiton scattering is considered. 
The absolute values of the spin correlation function between 
NN sites \{Si • Si+s)\ are 1 (bold lines), 0.5 (broken lines) and 
(dotted lines). Both the incident and scattered photon po- 
larizations are parallel to the z axis. Other parameter values 
are the same as those in Fig. 1. 
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with ut = (e 2 /mc 2 ) 2 and a bond length a. Pw is the 
polarization factor given by 



Pw = (e fei Ai); i.e kiXi ) v {e kf x f ) l [?k,\,) v 



(23) 



and S u (t) is the dynamical correlation function defined 

by 



S ll '(t) = (K l (t)lC 1 ' (Q)), 



with 



£'=£££«r(M + *0- 



(24) 



(25) 



TLr{i, i + Si) is a term of the effective Hamitonian given in 
Eqs. (3) and (15) concerning with a bond connecting site 
i and site i + Si. Here, l(= x, y, z) indicates a direction of 
the bond and 61 = ±1 The index r(= T u T 2 , A 1 ,E) clas- 
sifies the intermediate states. The exchange interaction 
Jr = t^/AEr in Hr is replaced by t%/(AEr — (w» — w/)) 
in TLr{i, i + Sf). This expression is obtained from the sec- 
ond order processes of the interactions between photons 
and electronic currents between the NN TM sites. 

We introduce the orbiton operators ip(k) with the ener- 
gies E(k) which are obtained by diagonalizing the Hamil- 
tonian represented by the Holstein-Primakoff boson op- 
erators tp(k). This Bogoliubov transformation is given 
by the matrix V(k) as 



Va(£) = £lMfe)VWfc). 



(26) 



ip(k) has 2N components in the system where the number 
of the OW mode is N, and ip a (k) (a > N) is the creation 
operator with the condition Va(^) = ^a-Ni—k)^ . For 
example, 

$(k) = { ai (k), ■ ■ ■ a N (k), a\(-k), ■ --aU-k)} ■ (27) 
By using the oprators, we obtain 

jc 1 = £ Mk)^h l aP (k)Mk) 



k,a,(3 



+ 



^£(^(0)+^>i(0)), (28) 



where h l a/3 (k) and g l a are the coefficients. In the pure 
00 states, the second term vanishes, i.e. the one-orbiton 
scattering is prohibited. Then the Fourier transform of 
the dynamical correlation function in the two-orbiton Ra- 
man scattering is given by 

2N N 

S»V) = AT£ £ £5{ W -£7 a (A)-£7 /3 (fc)} 

% a =N+ll3=l 

x \h l al3 {k)h l p a (k) + h l aj3 (k)h l a _ N p +N {%)). (29) 

We neglect the orbiton-orbiton interaction for simplic- 
ity, and S u (w) is represented by the convolution of the 
two OW modes. On the other hand, in the one-orbiton 
scattering, the dynamical correlation function reflects the 
OW at the momentum k = as 



JV 



S u '(u)=4Nj29U5W-E a (0)}. 



(30) 



In Fig. 4, the Raman spectra by OW in the (S-C/O- 
G) phase for LaV03, and in the (S-G/O-C) phase for 
YVO3 (the low temperature phase) are presented. The 
two-orbiton scattering processes are considered. In spite 
of this processes, a sharp peak structure appears at the 
lower edge of the continuum. This is attributed to the 
one-dimensional character of the OW and the factor 
h l a p(k) in Eq. (28) which enhances the lower edge. As 
for the selection rule, the Raman scattering is only ac- 
tive for the (zz) polarization where both the incident and 
scattered light polarizations, e*£ A and f % / ^ /J ar e paral- 
lel to the z axis. Through the interaction with the z- 
polarized lights, two electrons are exchanged between the 
d yz (d Z x) orbitals in the NN sites along the z axis. We 
mention that the local modes discussed in the previous 



section, i.e. the d x 



dy Z and d X y 



d zx excitations, 



are not detected by the Raman scattering, since the ex- 
change processes do not occur between the d xy and d yz 
(d zx ) orbitals. In Fig. 5, we show the Raman spectra 



G 




FIG. 5. Raman spectra 

for the {d y ( x+z )/ d y ( x - z ) / d x ( y+z )/ d x ( y - z ))-tyve 00 state with 
FM order for YTi03. The one-orbiton scattering is consid- 
ered. Symbols with vertical lines are the spectral weights of 
the delta- function. Continuous lines indicate S * f (uj) where 
the delta-functions are replaced by the Lorentz functions and 
the width of the function is chosen to be 0.25 Jo. The po- 
larizations of the incident and scattered photons are denoted 
as (A, A'). Other parameter values are the same as those in 
Fig. 2. 



from OW in the (d y{x+z) /d y ( x _ z) /d x{y+z) /d x(y _ z) )-typc 
00. The one-orbiton scattering is considered in the cal- 
culation. The spectra being active for the (xx) and (x'x') 
polarizations are the A\ g modes, and those for the (xx) 
and (x'y') ones are the B\ g ones. Here, the x, y and z 
directions are chosen to be parallel to the TM-0 bonds, 
and x' = -^( x + y) an d v' = 773 (~ x + 2/)- ^ is worth 
noting that all modes are inactive for the (zz) polariza- 
tion in contrast to the case of vanadates. This originates 
from a cancellation from the one-orbiton scattering con- 
tributions from site i and its NN site i + 5 Z along the 
z direction where the occupied orbitals have the mirror 
symmetry in terms of the xy plane between the two. 



B. inelastic neutron scattering 

Although the Raman scattering is a possible probe to 
detect OW as shown in the previous section, its observa- 
tion is limited to be the OW at zero momentum in the 
one-orbiton scattering, and the joint density of states of 
OW in the two-orbiton scattering. In the ti g 00 systems 
of our present interest, it is possible to detect by the 
inelastic neutron scattering. Formulate the differential 
scattering cross section in the scattering of initial (scat- 
tered) neutron with momentum fci(fc/), energy u>i(u)f), 
and polarization k(lf) (= x,y,z). The scattering cross 
section is given by 



"5 



(S-C/O-G) 



(S-G/O-C) 



'o.5 

n 



|<SS>|=0 



|<SS>|=1 




|<SS>|=1 



|<SS>|=0 



FIG. 6. Contour map of the scattering intensity S zz (ui) 
in the inelastic neutron scattering in the (S-C/O-G) and 
(S-G/0-C)-phases for LaV0 3 and YV0 3 (the low temper- 
ature phase), respectively. The intensity is plotted as a unit 
of 2/ Jo. The spin polarizations of both the incident and scat- 
tered neutrons are chosen to be parallel to the z axis. The 
reciprocal vector is G = (000). Other parameter values are 
the same as those in Fig. 1. 



da 2 



( 7 e2 



dflduif \rriNC 2 J \2' 



;gF(K) 



(31) 



with K — ki — kf, lj = u>i — ujf, and k = K/\K\. 
ghif ( w ) j g t nc F 0UT ier transform of the dynamical cor- 
relation function for the angular momentum operators 
L\ d defined by 

S l ^(f id -f i , d ,,t) = (L'- 1 (t)L l J d ,(0)), 

= 2{O idTl i i {t)O i , d , Tl i ! m, (32) 

where fa is the position of the d-th. TM ion in the i-the 
unit cell. In the linear spin wave approximation, S(K, ui) 
is obtained as 

N 



dd' a—1 
h* n l f JG-S dd , 



da a ™ 



(33) 



where D l da is defined by the Fourier transform of the 
angular momentum operator 



O dTll (k) = J2 D L(k)Mk). 



(34) 



S d d' is a vector connecting the d- and d'-th. TM ions in 
the same cell, and G is the reciprocal lattice vector. In 
contrast to the Raman scattering, the momentum de- 
pendence of the dispersion relation is detectable. The 
magnitude of the scattering intensity is expected to be 
the same order with that for the magnetic neutron scat- 
tering in magnets with L = 1. 
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°z r m xz r m x 
FIG. 7. Contour map of the inelastic neutron scattering 
spectra S" (ui) for the {d y{x+z) /d y{x _ z) /d x ( y+z) /d x{y _ z) )-type 
00 state with FM order for YTiOa. The intensity is plotted 
as a unit of 2/ Jo. The spin polarizations of the incident and 
scattered neutrons are chosen to be (a) unpolarized, and (b) 
parallel to the z axis. The reciprocal vector is G — (000). 
Other parameter values are the same with those in Fig. 2. 

We present the scattering intensities of the inelastic 
neutron scattering in the (S-C/O-G) and (S-G/O-C)- 
phases (Fig. 6), and those in the FM order with the 
{dy(x+z)/dy( x - z )/d x ( y+z )/d x ( y - z < ) )-type OO (Fig. 7). It 
is noted that the dispersive OW modes in vanadates are 
only detected by the z polarized neutron; the angular mo- 
mentum L z induces the excitations between the d yz and 
d zx orbitals. The local modes discussed in the previous 
section are active for the x and y polarized neutrons (not 
shown in Fig. 6). Therefore, the modes of OW are identi- 
fied by utilizing the polarized neutron scattering. This is 
also seen in the contour map of the scattering intensity in 
YTiOs; the scattering intensity for the OW modes in the 
higher energy bands are remarkable in the z-polarized 
neutrons. This implies that these modes mainly consist 
of the excitations between the d zx and d yz orbitals as 
explained in Sec. III. 

V. SUMMARY AND DISCUSSION 

In this paper, a theory of OW in perovskite titanates 
and vanadates with the triply degenerate ti g orbitals is 
present. We examine the dispersion relations of OW in 
the (S-C/O-G)- and (S-G/O-C) phases for LaV0 3 and 
YVO3 (the low temperature phase), respectively and 
that in the (d y{x+z) /d y{x _ z) /d x(y+z) / d x(y _ z) )-type 00 
for YTi03. We demonstrated that characteristics of OW 
in these systems can be detected by utilizing the Raman 
and inelastic neutron scatterings. 

In comparison with the OW in the manganites where 
the e g orbital is ordered, both the excitation spectra and 
the observation methods are distinct qualitatively in the 



present ti g orbital system. In particular, this is remark- 
ably seen in vanadates where the so-called pure orbitals 
are ordered. Thus, the selection rules for the Raman 
and neutron scatterings are strict. For example, the two- 
orbiton scatterings with the z polarized photons domi- 
nate the Raman spectra. This is attributed to the or- 
thogonality of the electron transfer integral between the 
NN vanadium sites. In the actual vanadates, there is the 
GdFe03-type lattice distortion which may make the one- 
orbiton scatterings possible. In the recent Raman scat- 
tering experiments in -RVO3, a new peak appears around 
60mcV in the (S-C/O-G) phase. 40 It is confirmed that 
this peak is active in the (zz) polarization configuration. 
We expect that this peak originates from OW excited by 
the two-orbiton scattering processes as discussed in Sec. 
Ill (Fig. 4). 

In the case of YTi03, the dispersion relation of OW 
and the Raman and neutron scattering spectra are more 
complicated than those in vanadates. This is because of 
the mixed 00 state with the four different orbitals in a 
unit cell. The present results are also distinct from those 
proposed in Ref. 18; the orbital excitations are exam- 
ined in the 00 states with high symmetry being differ- 
ent from the (d y ( x+z) /d y ( x ^ z) /d x ( y+z) /d x ( y ^ z) )-type 00 
and incompatible with the crystal symmetry of YTi03. 
In the present results, as shown in Fig. 2, it is found, 
in contrast to the previous results, 18 that there are the 
anisotropy of the dispersion relations in the xy plane and 
along the z axis, the two kind groups of the OW with 
higher and lower energies appear, and the flat bands are 
not seen along the (7r7T7r)-(7T7rO) direction. Actually, the 
inelastic neutron scattering experiments have started in 
YTi03. 41 ' 17 The detailed comparison between the theo- 
retical calculations and the experimental data can reveal 
nature of OW as well as that of 00 in the t 2g orbital 
systems. 
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